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Abstract 

Electron spin echo envelope modulation (ESEEM) has been observed for the first time from a 
coupled /leiero-spin pair of electron and nucleus in liquid solution. Previously, modulation effects in 
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, spin echo experiments have only been described in liquid solutions for a coupled pair of homonuclear 

o ■ 

spins in NMR or a pair of resonant electron spins in EPR. We observe low-frequency ESEEM (26 



and 52 kHz) due to a new mechanism present for any electron spin with S" > 1/2 that is hyperfine 
coupled to a nuclear spin. In our case these are electron spin (5 = 3/2) and nuclear spin (/ = 1) 



qh, in the endohedral fullerene NQCgo- The modulation is shown to arise from second order effects in 



the isotropic hyperfine coupling of an electron and ^^N nucleus. 
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I. INTRODUCTION 



Measuring the modulation of a spin echo in pulsed magnetic resonance experiments has 
become a popular technique for studying weak spin-spin couplings. It is used extensively in 
the fields of chemistry, biochemistry, and materials science, both in liquids and sohds, using 
nuclear magnetic resonance (NMR) and electron paramagnetic resonance (EPR) ^ O, . 
Two distinct mechanisms for spin echo modulation have been identified in the literature. 

In the first mechanism, a pair of spins, S and J, are coupled through an exchange or 
dipole-dipole interaction, J-S-I or J-Sziz, and the echo modulation arises for non-selective 
refocusing pulses which flip both coupled spins. The magnitude of the echo signal oscillates 
as cos(Jt), where t is the interpulse delay time page 500). This is most commonly 
observed for coupled pairs of homonuclear spins 5], though it is also known for pairs of 
coupled electron spins 

It should be emphasised that a similar coupling, J-S-I, between unlike spins (including 
coupling between heteronuclear spins, hyperfine coupling between electron and nuclear spins, 
and electron-electron coupling between two electron spins with different Larmor frequencies) 
results in no modulation effects from this mechanism: The hetero-spin coupling energy 
changes its sign upon application of the refocusing pulse (because only one spin flips), and 
the magnetization is thus fully refocused at the time of echo formation, in the same way as 
in the presence of any other inhomogeneous magnetic fields. 

A second ESEEM mechanism, which does apply to coupled pairs of hetero-spins, has 
also been identified . This mechanism requires anisotropic spin-spin interactions (e.g. 
^zzSzh + ^zxSzIx) and is therefore restricted to solids or high viscosity liquids. The modu- 
lation arises as a result of "branching" of the spin transitions created by the refocusing pulse. 
The resonant spin S precesses with Larmor frequency that is different before and after the 
refocusing pulse and therefore accumulates an additional phase which causes oscillations in 
the echo signal, as cos{uJikt), where uik is the spin transition frequency of the non-resonant 
spin I. The amplitude of the oscillations depends on magnitude of the anisotropic hyperfine 
component. 

In this Paper we demonstrate that, contrary to previous belief, echo modulation effects 
can also be observed for a hetero-spin pair coupled by a purely isotropic spin interaction, and 
we thus identify a new ESEEM mechanism. Our hetero-spin pair is the endohedral fuUerene 
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NOCeo in CS2 solution, with electron spin S = 3/2 interacting through an isotropic hyperfine 
coupling {a-S-I, a = 15.8 MHz) to the nuclear spin / = 1 of ^^N. The isotropic hyperfine cou- 
pling lifts the degeneracy of the electron spin transitions, leading to a profound modulation 
of the echo intensity at about 52 kHz. We shall show that this third modulation mecha- 
nism is only effective in high-spin electron systems {S > 1/2). The N^Cqo molecule has an 
exceptionally long electron spin dephasing time (T2 = 210 /xs), enabling the observation of 
this low frequency ESEEM for the first time. 

II. MATERIALS AND METHODS 

High-purity endohedral NtSCeo was prepared 0], dissolved in CS2 to a final concentration 
of 1-2-10^^/cm^, freeze-pumped in three cycles to remove oxygen, and finally sealed in a 
quartz EPR tube. Samples were 0.7-1.4 cm long, and contained approximately 5 ■ 10^^ 
NQCgo spins. Pulsed EPR measurements were done at 190 K using an X-band Bruker 
Elexsys580e spectrometer, equipped with a nitrogen-flow cryostat. In the 2-pulse (Hahn) 
electron spin echo (ESE) experiments, 7i/2 — t — ti — t — echo, the tt/2 and vr pulse durations 
were 56 and 112 ns respectively. Phase cycling was used to eliminate the contribution of 
unwanted free induction decay (FID) signals. 

III. RESULTS AND DISCUSSION 

Fig. CJA) shows the continuous-wave EPR spectrum of N@C6o in CS2 at room tempera- 
ture. The spectrum is centered on the electron g-factor g = 2.0036 and comprises three lines 
resulting from the hyperflne coupling to ^^N Q]. The relevant isotropic spin Hamiltonian (in 
angular frequency units) is: 

Ho = UeS^ - uih + a-S-I, (1) 

where Ue = gPBo/h and uj = gj^nBo/h are the electron and ^^N nuclear Zeeman frequencies, 
g and gi are the electron and nuclear g-factors, (3 and /?„ are the Bohr and nuclear magnetons, 
h is Planck's constant and Bq is the magnetic field applied along z-axis in the laboratory 
frame. Each hyperfine line (marked in Fig. ^A) with M/ = and ±1) involves the three 
allowed electron spin transitions AM^ = 1 within the 5 = 3/2 multiplet. These electron 
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FIG. 1: (A) EPR spectrum of NOCgo in CS2 at room temperature. Each line in the triplet signal 
is labeled with the corresponding projection Mj of the ^^N nuclear spin. (B-D) Zoom-in for each 
line showing details of the lineshape structure. Small satellite lines (marked with *) are due to 
a hyperfine interaction with the natural abundance of ^'^C nuclei on the Ceo cage. Measurement 
parameters: microwave frequency, 9.67 GHz; microwave power, 0.5 ^W; modulation amplitude, 
2 mG; modulation frequency, 1.6 kHz. 

spin transitions remain degenerate for Mj = as seen in Fig. ^C) but split into three lines 
(with relative intensities 3:4:3) for Mj = ±1, as seen in Figs. ^B) and (D). This additional 
splitting of 0.9 fiT originates from the second order hyperfine corrections a? /oJe = 26 kHz, 
and its observation is only possible because of the extremely narrow EPR linewidth < 0.3 
in NtQCgo- Similar second-order splittings have been reported for the related spin system of 
endohedral fuUerene ^^P@C6o which has S = 3/2 coupled with / = 1/2 [lol |. 

Fig. 12^ A) shows two-pulse echo decays measured at the central Mj = and the high-field 
Mj = —1 hyperfine lines. The decay is monotonic for Mj = and has an exponential 
dependence exp(— 2r/T2) with T2 = 210 ^s. However, the decay is oscillatory for Mj = —1 
(and also for Mj = +1, not shown) — the Fourier transform of the decay reveals two peaks 
at frequencies 26 and 52 kHz as seen in Fig. EfB). These frequencies correlate closely to the 
splitting of 26 kHz found in the EPR spectrum in Fig. E^B) and (D), indicating that the 
two effects have the same origin. 

We shall use the spin density operator formalism to derive the modulation effects for the 
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FIG. 2: (A) Two-pulse ESE decays for NQCeo in CS2 at 190 K measured at the central Mj = and 
the high- field Mj = — 1 hyperfine components of the EPR spectrum. (B) The Fourier Transform 
(FT) of the oscillatory echo decay at Mj = — 1. 



spin system S = 3/2, I = 1. The spin density matrix after our two-pulse echo experiment 
is given by 



Ct{t) = {UrR2UrRl) " (Tq " {UrRlUrRl)^ . (2) 

Here, (Tq is the density matrix at thermal equilibrium; in the high-temperature approx- 
imation valid in our experiments, ao can be substituted with a spin operator 5*^ j^. The 
evolution operator Ur = exp(— zTYot) describes a free evolution of the spin system between 
the applied microwave pulses, and the spin rotation operators, R^, describe spin rotation 
upon application of the two microwave pulses, i = 1,2. The measured echo intensity is 

Vir) = Tr [a(r) ■ D] . (3) 

The detection operator D = Sy ^ Vmj involves the ^^N nuclear spin projection operator 
Vmi to selectively detect only those spin transitions associated with a specific nuclear spin 
projection M/. In a pulsed EPR experiment, this corresponds to performing measurements 
at the resolved hyperfine line in the EPR spectrum of NQCeo and integrating over the echo 
signal shape to average out oscillating signals from other off-resonance hyperfine lines. The 
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echo derives from a sum of the single quantum (SQ) coherences (represented by the terms 
o"n.,n+i and o"„+i,„ in the density matrix), weighted by factors from the detection operator D. 

Our Hamiltonian, TYq (Eq- D), has small off-diagonal elements provided by the a{IxSx + 
lySy) terms. These terms are often omitted since they are known to contribute only sec- 
ond order energy corrections P, 11|, but are, in fact, directly responsible for the observed 
ESEEM. Diagonalisation of Eq. yields the magnitude of these corrections to be of or- 
der 5 = a^/uJei consistent with the splitting observed in Fig. ^B) and (D). In addition, 
the isotropic hyperfine interaction introduces a small degree of mixing between J^, Sz basis 
states, however we find that this small mixing need not be considered to appreciate the ori- 
gin of the observed ESEEM. With this assumption, our Hamiltonian and all other operators 
have block- diagonal structures with non-zero elements only between states with the same 
Mi. Transitions with simultaneous flip of both electron and nuclear spins are thus forbid- 
den and the evolution of electron spin can be treated individually for each nuclear spin 
manifold. We can therefore avoid the derivation in the full 12 x 12 Hilbert space in a general 
form, and instead reduce the dimensionality to 4 x 4. The validity of this approximation is 
confirmed below by a rigorous derivation using Average Hamiltonian Theory. The reduced 
Mi = +1 subspace of the diagonalised Hamiltonian, correct to second order in a, becomes: 
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The term —Iz{uji + 7/4:6) represents a constant energy shift and can be ignored. We move 
into a resonant rotating frame (the coordinate system rotating with the microwave frequency 
ujjnw around the laboratory 2;-axis). In this frame the spin Hamiltonian (j3]) transforms to 
'H — UmwSz, such that: 



6 



\ 



(6) 



3/2A + 5 

1/2A 
-1/2A 

y -3/2A + S J 

where A = uJe + a + 6/2 — is the resonance offset frequency. 

The rotation operator = exp{—i(Ho + Hi)tpJ can be simphfied by taking into account 
a finite excitation bandwidth of the microwave pulses, which are selective upon one hyperfine 
line in the EPR spectrum. Tii = gjSBiSx/h where Bi is the microwave magnetic field applied 
along X-axis in the rotating frame and tp. is the duration of the microwave pulses. Further- 
more, since the bandwidth of Bi ^ 4.5 MHz is large compared to both the intrinsic EPR 
linewidth 9 kHz and the second-order spitting 26 kHz of the outer lines in the EPR spec- 
trum, all three spin transitions within electron the S = 3/2 multiplet are equally excited and 
the respective rotating operator can be approximated as ~ exp(— zTiitp.) = exp{—i9iSx), 
where 6i = g(3Bitp./h is the microwave pulse rotation angle. This results in the following 
spin rotation operator for an on-resonance hyperfine line: 
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We need not consider the rotation operator for off-resonance lines, as their excitation will 
lead to an oscillating echo signal which is averaged out by selective detection (i.e. through 
the detection operator D) and therefore does not contribute to the overall echo signal. 

We evaluate Eq. |21for the two-pulse echo experiment 7r/2 — r — vr — r — ec/io, and rearrange 
terms for the purposes of the discussion which follows: 
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Evaluating the heart of the echo sequence, UrR^Ur, is instructive in understanding the 
source of the observed modulation. From Eq. [71 a perfect vr rotation is: 
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with the resulting echo sequence operator shown below. 
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This indicates that over the course of the experiment, states S = ±3/2 pick up a phase 
of 26t with respect to the states S = ±1/2, or in other words, the outer SQ coherences 
{<^i,2, cr2,i, cr3,4 and (74^3) oscillate with frequency 26, while the phases of the inner SQ coher- 
ences (cT2,3 and (73^2) remain constant. 

The initial SQ coherences are provided by the first rotation (vr/2): 



Rn/2 ■ ^0 ■ {Rl/2)'^ 
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(11) 



i 

y iy3/2 

and the measuring weighting factors for each coherence are: V3/2, 1, and ^3/2 as- 
sociated with the detection operator V). Together these imply that the three coherences 
contribute to the measured echo intensity with relative amplitudes 3:4:3. In other words, 
their sum will yield a constant component, and one oscillating with frequency 25, with 
respective amplitudes 4:6. This is confirmed upon evaluation of Eq. IHl 



Vmj=±i{t) =2±3cos25r. 



(12) 



and is consistent with the observed echo in Fig. EB- The modulation amplitude is deep and 
the echo signal can change its sign at the minima. 

The effect is illustrated in Fig. E^A), which shows the phases gained during the "defo- 
cusing" period r, i.e. free evolution after the initial pulse Oi and before the refocusing pulse 
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FIG. 3: Phases gained by the density matrix elements during the 2-pulse ESE experiment at 
Mj = +1. Phases caused by off-resonance refocus fully and are therefore omitted for clarity. 
The shaded off-diagonal elements represent single-quantum (SQ) coherences which generate the 
echo signal upon refocusing. In (A) the arrows indicate the transition between the SQ coherences 
caused by a perfect refocusing pulse with 62 = vr. (B) shows all possible transitions (spin coherence 
branching), for one SQ coherence element, caused by an imperfect refocusing pulse with 62 7^ it. 



62- These phases derive from the differences between adjacent elements along the diagonal 
of the Hamiltonian in Eq. El (the off- resonance. A, is ignored as it is fully canceled upon 
echo formation). The six SQ coherence elements, shaded in the Figure and responsible for 
echo formation, gain the phases or ±6t. Upon application of the perfect refocusing pulse 
with 62 = 7T each SQ coherence element uniquely transforms, as shown with arrows, and 
continues to evolve during the "refocusing" period r to gain an additional phase which does 
not compensate, but instead doubles, the initial phase. Thus, at time of echo formation the 
SQ coherences arrive with three different phases and ±26t. Their vector sum interferes 
destructively to produce an echo signal whose magnitude oscillates as 26t in accordance 
with Eq. d 

The case of an imperfect refocusing pulse with 62 ^ it is shown in Fig. El^B). In contrast 
to 62 = TT with a one-to-one transformation of each density matrix element, the non-ideal 
pulse generates branching of the electron spin transitions. Therefore, the SQ coherence 
element which initially gained the phase +6t during the "defocusing" period, refocuses into 
three SQ coherences (shown with the arrows), each accumulating different phases during the 
"refocusing" period. Thus, at the time of echo formation the accumulated phases are 0, +6t, 
and +26t. The vector sum of these and other SQ coherences produces a complex interference 
with the echo signal oscillating with two frequencies 6 and 26 as observed in Fig. |2l Thus, 
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the second harmonic 5 found in the echo modulation is the result of an imperfect refocusing 
pulse. 

The preceeding physical description provides an intuitive view of the new ESEEM effect. 
For it to be fully rigorous, however, we should consider the effect of the mixing of the J^, 
basis states. This necessarily involves the full 12-dimensional Hilbert space and the argument 
rapidly becomes opaque. Fortunately, the same results can be rigorously obtained from the 
first order correction in Average Hamiltonian Theory (AHT) 13|, which is equivalent to a 
standard perturbation theory approach in the rotating frame. 

We begin by transforming the original Hamiltonian (Eq. (0)) into a rotating frame of 
angular frequency Umw, defining f2e as the deviation from the electron Larmor frequency, 

fig ^mwi 

Ho = ^eSz - i^llz + a[SzIz + Sxlx COs{uJmwt) 

-SyIySm{uJmwt)]. (13) 
The oscillatory terms in Eq. El are averaged out in the zeroth-order average Hamiltonian: 

n'o^ =n,Sz-uJiIz + aSzIz. (14) 

Here, the bar over Hq refers to an average over one period of the oscillation Umw As the 
Hamiltonian in Eq. (jl4p results in no modulation of the echo signal for a coupled hetero-spin 
pair, the higher order terms of the average Hamiltonian must be included (see, for example, 
Ref. |3|, p. 83). The first order correction is: 



/(/ + 1) - Sz - (SiS + 1) - S!) Iz . (15) 



;t7(1) ^ 
=2 

We find that the average Hamiltonian Tio = '^o'^^ + "H^o^ is sufficient to describe the modula- 
tion effects. In this approach the time-dependent mixing terms have been properly accounted 
for to produce (after time averaging) the second order energy corrections in Tig^''. However, 
these mixing terms appear to average to zero, with the result that the effective Hamiltonian, 
Tio, is a diagonal matrix, thus validating our earlier qualitative approach. 

It can also be verified that in the presence of the applied microwave field, the average 
Hamiltonian is, to first order, the simple sum Tig + gPBiS^/h. Therefore, the rotation 
operator, Rf, and detection operator, D have the same block diagonal structures described 
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above (non-zero elements only within the same Mj subspace). This allows us again to reduce 
the dimensionality of the Hilbert space to 4 x 4. 

Substituting Eqs. (I7[ HH and HH)) into Eq. (jH)), and after some manipulation, we find the 
following expressions for the echo amplitude in a general two-pulse sequence 9i — t — 62 — 
T — echo, for the S = 3/2, 1 = 1 spin system. The echo modulation is identical for the two 
outer Ml = ±1 hyperfine lines, 

VM,=±i{r) = 2sin^isin2^[Ao(^2) + ^i(^2)cos5r 

+A2{92)cos26t], (16) 



where 



^0(^2) = l-6cos^| + ^cos^|, 

A,ie2) = 6 cos^ I (^2 - 3 cos^ 1^ , (17) 
MO.) = ^sin^|(l-3cos^|). 

The signal is modulated with frequencies S and 2S, consistent with the experimental obser- 
vations in Fig. 121 The modulation amplitudes in Eq. |T7jl depend strongly on the rotation 
angles of the microwave pulses. At optimal rotation angles, 9i = 71/2,62 = tt, Eq. IT^ is 
recovered. 

In contrast to the two outer lines, the echo signal at the central M/ = hyperfine line 
shows no modulation effects, 

FM,=o(r) = 2sin^isin^^. (18) 

It is instructive to consider which terms in the average Hamiltonian Tig^"* give rise to the 
modulation effects. The terms I{I + l)Sz and S{S + l)Iz are not responsible as they produce 
only a constant shift to the electron and nuclear Zeeman frequencies, respectively. The term 

stays invariant during the experiment. This term changes its sign, but not its magnitude, 
during the refocusing pulse and thus fully refocuses. Therefore, SlI^ is solely responsible 
for the modulation effects. For M/ = 0, this last term becomes zero and consequently there 
are no modulation effects produced at the central hyperfine line in the EPR spectrum. Note 



I^Sz is also irrelevant because electron-nuclear flip-fiop transitions are forbidden, hence 
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FIG. 4: Two-pulse ESE decays (A) and their Fourier Transform spectra (B) measured at hyperfine 
line Mi = — 1 of the EPR spectrum, using the refocusing pulse 02 = tt and 62 = 2tt/3. Other 
experimental conditions are the same as in Fig. [^J 

that this term is responsible for the spin- dependent shifts to the energies of the spin states 
used in the earlier derivation (see Eqs. El and EI). 

The effect of a non-7r refocusing pulse angle is shown in Fig. |3J As predicted by Eq. (fTB|) . 
we observe that when 62 = 27r/3, the modulation effects are dominated by the low-frequency 
6, rather than the high-frequency 26 found when 62 = tc. 

Eq.Elalso confirms that a perfect vr refocusing pulse yields only a 26 frequency component 
in the modulation, however, a 6 component is clearly observed in Fig. |21 The imperfection 
is explained by the inhomogeneity of the microwave magnetic field Bi in the resonator 
cavity which results in a distribution of spin rotation angles 62 across the ensemble. If we 
assume a Gaussian distribution of rotation angles, the relative intensities of the low- and 
high-frequency components /1//2 = 0.17 in the experimental spectrum corresponds to a 
standard deviation of cr = 0.31 radians. This corresponds to a 10% error in a vr rotation 
angle, consistent with previously reported values for i?i-field inhomogeneity in this resonator 
cavity Q]. To verify that the Bi field inhomogeneity is the source of the low- frequency 
component, we applied an error- correcting composite vr-pulse as the refocusing pulse. The 
resulting ESEEM contained only the single 52 kHz frequency component. 
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The derivation of the modulation is easily generalised to the case of an arbitrary electron 
spin S > 1/2 coupled through an isotropic hyperfine interaction to a magnetic nucleus. 
Using an approach similar to that described in Eqs. (jUl IIUI and lll|). the general expression 
for 2-pulse ESEEM with a perfect refocusing pulse can be shown to be 

s 

V{t) = (S- Ms) {S + Ms + 1) e'(i+2^^s)M,5r ^^g^ 

Ms=-S 

The summation is over electron spin projections Ms, whilst the nuclear spin projection 
Ml identifies the hyperfine line of the EPR spectrum in which the modulation effects are 
observed. 

IV. CONCLUSIONS 

Potential applications of this new mechanism include measuring the hyperfine coupling 
constant and determining electron spin number; it may also be relevant to certain quantum 
information processing schemes . The accurate measurement of the hyperfine constant in 
a continuous- wave EPR measurement is subject to Bq field instability (typically >10 mG). 
However, the ESEEM frequency can be measured accurately, given a sufficiently long de- 
phasing time, potentially providing a more precise measurement. In this case accuracy may 
be improved by moving to lower applied magnetic fields (lower EPR frequency). In con- 
trast with other types of ESEEM, this would also lead to a higher modulation frequency. 
Finally, we note that given appropriate electron-nuclear spin coupling energies and decoher- 
ence times, this effect will also lead to modulation in nuclear spin echo experiments. 
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